Spin and electronic correlations in gated graphene quantum rings 
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We present a theory of graphene quantum rings designed to produce degenerate shells of single 
particle states close to the Fermi level. We show that populating these shells with carriers using 
a gate leads to correlated ground states with finite total electronic spin. Using a combination of 
tight-binding and configuration interaction methods we predict ground state and total spin of the 
system as a function of the filling of the shell. We show that for smaller quantum rings, the spin 
polarization of the ground state at half filling depends strongly on the size of the system, but reaches 
a maximum value after reaching a critical size. 



I. INTRODUCTION 

There is currently significant interest in developing un- 
derstanding of electronic properties [l|-y, lD, |8y and ap- 
plications [5|-[7| of graphene. Starting with graphene as 
a zero gap non-magnetic material, reducing the lateral 
size and controlHng the shape and character of the edges 
opens the possibility of controlling the energy spectrum 
and hence electronic and magnetic properties of graphene 
nanostructures [1Q|414J. In particular, the zigzag edges 
are responsible for degenerate energy shells at the Fermi 
level [13|, I15M19I I22I , and associated finite spin polariza- 
tion as a result of electron-electron exchange interactions 
I5I423J. However, the coupling of spin polarized zigzag 
edges was shown to be anti-ferromagnetic in graphene 
nanoribbons [15] with no net spin polarization. In con- 
trast, in triangular graphene quantum dots spin polarized 
edges were shown to couple ferromagnetic ally leading to 
a finite magnetic moment |19l42ll \. 



. The purpose of 
this work is to answer the question whether it is possible 
to use graphene nanoribbons to build graphene nanos- 
tructures with finite magnetic moment. We show here 
that by designing a hexagonal ring from six ribbons we 
obtain a quantum system with degenerate shells in the 
energy spectrum. By filling these shells with additional 
electrons using metallic gate we obtain maximally spin 
polarized ground state for the half filling of the degener- 
ate shell. 

Semiconductor quantum rings have been investigated 
by a number of groups 2j427|. The ring geometry al- 
lows to observe quantum phenomena, e.g. persistent cur- 
rent and quantum interference effects |28|], in particular 
Aharonov-Bohm (AB) oscillations [2^. The AB oscilla- 
tions manifest themselves as periodic oscillations in the 
energy spectrum of the electronic system as a function 
of the number of flux quanta entering the ring [25]. The 
AB effect for a single electron in single lithographically 
defined semiconductor quantum ring [26L hole in a type- 
II semiconductor dot [27] and exciton [29|, |30| in a finite 
ring, was demonstrated. 

The Aharonov-Bohm conductance oscillations were 



also recently observed in a graphene ring [3l|, |32|. 
The electronic properties of a single Dirac Fermion in 
graphene quantum rings were studied using effective mass 
[33, 34] and tight-binding methods [33, 35, 36]. Valley 
degeneracy in graphene was shown to be lifted by the 
magnetic field [33, 35] since the magnetic field has the 
opposite sign in the two valleys. 



In this work we combine the tight binding method with 
the configuration method to determine the electronic and 
spin properties of graphene quantum rings with zig-zag 
edges as a function of the number of additional electrons 
controlled by the gate. We analyze the energy spectrum 
of quantum rings as a function of the width of the ring 
and its size. We find degenerate electronic shells near 
Fermi energy for the thinnest structures. We use config- 
uration interaction method to treat exactly interaction of 
additional electrons in the degenerate shell as a function 
of shell filling. We determine the dependence of the spin 
polarization of the ground state as a function of the fill- 
ing of the degenerate shell and the size of the structure. 
We show that by changing the size of the structures we 
control the splitting between levels in a degenerate shell 
which in turn significantly influences magnetic properties 
of the ground state. The stabilization of the spin phase 
diagram at a critical size of the ring is observed. 



This paper is organized as follows. In Section II, we 
introduce both the tight-binding (TB) model for single 
particle levels and configuration interaction method for 
electron-electron interaction. In Section III, we show the 
method for constructing hexagonal ring structures with 
different width and length. Next, in section IV we present 
discussion of single particle energy spectra and explain 
the origin of the shell structure. In section V we ana- 
lyze the effect of electron-electron interactions and spin 
properties of electrons in degenerate shells as a function 
of shell filling. Finally, in section VI we summarize ob- 
tained results. 



II. MODEL AND METHOD 

The single particle energy spectrum of Hz electrons 
in graphene quantum rings can be obtained using the 
tight-binding Hamiltonian [37]. The tight-binding model 
was successfully applied to carbon materials such as 
graphite, graphene, nanoribbons, nanotubes, fullerens, 
and graphene quantum dots [13, 14, ll6Ml9L l37| . The 
Hamiltonian in the nearest neighbors approximation can 
be written as 



H=t^ 



cLc 



a^jcr^ 



(1) 



{hJ),o- 



where t is the hopping integral, q^ and Cia- are creation 
and annihilation operators of electron on II^ orbital on 
site i with spin a =t,i respectively, and (z,j) indicate 
summation over nearest-neighbours. Diagonalization of 
the tight-binding Hamiltonian generates single particle 
energies e^ and single particle orbit als |s,cr). 

In order to include electron-electron interactions the 
many-body Hamiltonian Hmb is written as 

+ 2 5Z ^^^ I ^ I df)al^al^,ada'afa, (2) 

where the first term describes single particle energies ob- 
tained from the tight binding Hamiltonian given by Eq. 
(dl and the second term describes interactions between 
particles occupying these single particle states. By using 
Slater tt^ orbitals [38] we calculated two-body Coulomb 
matrix elements (zj \V \kl)^ where i^j^k^l are the site 
indices. In numerical calculations, on-site and all scat- 
tering and exchange terms up to next nearest neighbor's 
are included. Few largest Coulomb matrix elements are 
given in the Appendix. We use t = — 2.5eV for the hop- 
ping integral, and the effective dielectric constant is set 
to K, = 6. 

III. CONSTRUCTION OF GRAPHENE RINGS 

In this section we show the method for constructing 
a mesoscopic ring from six graphene nanoribbons. Such 
a choice of building blocks helps us to understand single 
particle spectrum of the ring. In Fig. [T]we show two sets 
of six graphene ribbons arranged in a hexagonal ring. 
Each ribbon consists of two types of atoms from the unit 
cell of honeycomb lattice, indicated by red (light grey) 
and blue (dark grey) circles in Fig. [H On the left side, 
thinnest possible ribbons with one benzene ring width 
are shown, denoted as W^ = 1. Each of them consists of 
16 atoms (the length L = 4, is measured by the number 
of one type of atoms in the upper row), so the final ring 
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FIG. 1. (Color online) The method of construction of ring 
structures from smaller ribbon-like units. On the left, there 
are six thinnest possible ribbons (one benzene ring thick, what 
we note as VK = 1) arranged in a hexagonal ring structure. 
The length of each ribbon is given by L = 4, the number of 
one type of atoms in one row. Each ribbon consists of 16 
atoms which gives a total of 96 atoms in a ring. On the right, 
there are six ribbons with width W — 2 (two benzene ring 
thick). Each of them consists of 21 atoms giving a total of 
126 atoms in a ring. We create a thicker ring with similar 
length (L = 4) but smaller antidot inside. 



is built of 96 atoms. Small black arrows indicate bonds 
(and hopping integrals between nearest neighbors in a 
tight-binding model) between neighboring ribbons (two 
arrows in the case of thinnest structures). The number of 
such connecting atoms increases with increasing width as 
seen on the right hand side of Fig. [H The thicker ribbon 
{W = 2) has identical length to the one from the left side 
(L = 4). In this case there are three connecting atoms 
(three small black arrows indicate three bonds). The final 
ring is built of 126 atoms. By connecting neighboring 
ribbons with different lengths and widths we create rings 
with different single particle spectra. 



IV. SINGLE PARTICLE SPECTRA 

In Fig. [2] we show the single-particle energy levels 
near Fermi level obtained by diagonalizing tight-binding 
Hamiltonian, Eq. ([T]), for rings with a given length 
(L = 8) and different widths. The thinnest ring (W = 1) 
consists of 192 atoms. For this structure we observe 
nearly degenerate shells of energy levels separated by 
gaps. Each shell consists of six levels: two single and 
two doubly degenerate states. First shell over the Fermi 
level is almost completely degenerate while in the second 
one the degeneracy is slightly removed. We note that for 
rings with different lengths, the gap between first and sec- 
ond shell is always larger then the gap at the Fermi level. 
With increasing width of the ring, the spectrum changes 
completely. For the rings with width W = 2 (270 atoms), 
W = 3 (336 atoms) and W = 5 (432 atoms) shells are 
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FIG. 2. (Color online) Single particle spectrum near Fermi 
level for ring structures with L = 8 and different widths. The 
shell structure is clearly observed only for the thinnest ring 
(W = 1). Dotted blue (grey) line indicate the location of 
Fermi energy. 



not visible. For W = 4 (390 atoms) we observe appear- 
ance of shells separated by gaps further from Fermi level 
but the splitting between levels in these shells is much 
stronger in comparison to the thinnest ring. We note 
that for 1^ > 2, although we do not observe a clear pat- 
tern of shells around the Fermi level, single shells of six 
levels separated by gaps from the rest of the spectrum 
appear far away from the Fermi energy in some cases. 

In order to have a better understanding of the struc- 
ture of the tight-binding spectra, in Fig. [3] we show the 
evolution of single particle energies from six independent 
ribbons to a ring as the hopping between the ribbons 
is increased. To achieve this, first we diagonalize tight- 
binding Hamiltonian matrix for a single ribbon. We then 
take six such ribbons and create Hamiltonian matrix in 
the basis of the eigenvectors of six ribbons. Here, the 
matrix has diagonal form. All energy levels are at least 
six fold degenerate. Next, using the six ribbons basis, 
we write hopping integrals corresponding to connecting 
atoms between neighboring ribbons (indicated by small 
black arrows in Fig. [T]). By slowly turning on the hop- 
ping integrals and diagonalizing the Hamiltonian at every 
step, we can observe the evolution of the spectrum from 
single particle states of six independent ribbons to a ring. 

The hopping integrals between connecting atoms of 
neighboring ribbons are indicated by f in Fig. [3l For 
the thinnest ring (Fig. [3fa)), each ribbon consists of 32 
atoms. There are only two connecting atoms between 
neighboring ribbons, giving only two hopping integrals t' 
between each two ribbons in the nearest neighbors ap- 
proximation. We see that their influence is very small 
and six fold degenerate states evolve into shells with a 
very small splitting between levels. We note that this 
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FIG. 3. (Color online) The evolution of the single parti- 
cle spectrum from six independent ribbons with L = 8 to a 
hexagonal ring structure spectrum, t' indicate hopping inte- 
grals between neighboring ribbons, a) For the thinnest ring 
(VK = 1) six fold degeneracy is slightly removed, preserving 
a shell structure. For thicker structures ((b) and (c), W = 2 
and W = 3 respectively) the six fold degeneracy is strongly 
lifted and shell structure is not observed. 



splitting is a bit stronger for higher energy levels but 
due to large gaps between consecutive levels of single rib- 
bon the shell structure is still clearly observed. For the 
thicker structures (Fig. ^h) and (c)) the evolution of the 
spectrum has a more complicated behavior. For a given 
ring each ribbon consists of different number of two types 
of atoms which give rise to zero energy edge-states [39| . 
With increasing width, the number of zero-energy states 
increases as well as the number of connecting atoms (and 
equally the number oft' hopping integrals). This causes a 
stronger splitting of levels for thicker rings in comparison 
to the thinnest one. Thus, the thicker ring's spectrum 
close to the Fermi level is due to the splitting of zero- 
energy states of independent ribbons. For W = 2 (one 
zero-energy state) and W = 3 (two zero-energy states), 
each ribbon consists of 45 and 56 atoms respectively, and 
the evolution of their spectrum is similar. The degener- 
acy is strongly lifted and no shell structure is observed. 
In order to illuminate the influence of f hopping in- 
tegrals on the thinnest ring spectrum, in Fig. |4]we also 
show the electronic densities for the first shell over the 
Fermi level for three different values of f (indicated in 
Fig. [31(a)). For t' = there are six independent rib- 
bons and first shell is perfectly six fold degenerate. The 
electronic charge density in each ribbon is larger on the 
two atoms with only one bond (see FigH]) and gradually 
decreases along the length. For t' = 0.5t the total en- 
ergy of the shell increases and the degeneracy is slightly 
removed. Here, the highest peak of electronic charge den- 
sity is moved towards the center of each ribbon in com- 
parison to t' = case. Increasing t' to t causes increasing 
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FIG. 4. (Color online) Energy levels and corresponding total 
electronic densities for the first six states over the Fermi level 
for the thinnest structure with L = 8 (192 atoms), for i) 
t^ = 0, ii) t^ = 0.5t, iii) t' = t. The three values of t' hopping 
integrals are indicated in Fig. (Sja). 
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FIG. 5. (Color online) Single particle spectrum near Fermi 
level for the thinnest ring structures with length L — 8 and 
L = 4. The shell structure is clearly observed. The splitting 
between levels in the first shell is smaller for larger structure. 
Dotted blue (grey) line indicate the location of Fermi energy. 



of the total energy of the shell and the highest peak of 
electronic charge density is now perfectly in the middle 
of each arm of the ring. Thus, both electronic charge 
density and energy of levels change slightly during the 
gradual transition of ribbons into a hexagonal ring struc- 
ture. 

We find degenerate shells near the Fermi energy only 
for the thinnest rings. In the rest of the paper, we will 
focus on the single and many particle properties of these 
structures as a function of their lengths. In Fig. [5] we 
show the low energy spectrum for two thinnest rings 
with different lengths. We clearly see shells with six lev- 
els. The splitting of levels from the first shell over the 
Fermi level is smaller for larger ring. For ring structure 
with L = 4 the difference between the highest and the 
lowest energy of levels forming the first shell is around 
0.069t ^ 0.17 eV. In comparison, for ring with L = 8 this 
value is around 0.006t ^ 0.015 eV. Thus, we conclude 



that for smaller rings single particle energies can play 
important role in the properties of many particle states 
while for the larger rings interactions are expected to be 
more important. 



V. ELECTRONIC INTERACTIONS IN A 

DEGENERATE ELECTRONIC SHELL OF A 

GRAPHENE QUANTUM RING 

In this section, we study ground and excited states as 
a function of the number of additional interacting elec- 
trons in degenerate shells of quantum rings with differ- 
ent size L. In our calculations we assume that all states 
below the Fermi level remain fully occupied. This is jus- 
tified as long as there is a sufficiently large energy gap 
at the Fermi level. Next, we add extra electrons to the 
charge neutral system. In a first approximation we ne- 
glect scatterings from/to the states below Fermi energy. 
Moreover, because of the large energy gap between first 
and second shell we can neglect scatterings to the higher 
energy states. This allows us to treat first shell of the 
thinnest ring as an independent system which signifi- 
cantly reduces the dimension of the Hilbert space. All 
the shells in the studied structures consist of six levels. 
For a given number of extra electrons we create a ba- 
sis of all possible configurations of electrons distributed 
within these states. Since the total spin of the system 
is conserved, we diagonalize the many-body Hamiltonian 
in subspaces with total projection of spin onto z-axis Sz- 
The largest dimension of the Hilbert space is for the half 
filling (6 electrons) and 5'^ = for which there are 400 
configurations. 

First, we study the magnetic properties of the half 
filled shell. Figure [6] shows the low energy spectra for 
the different total spin S of half filled first shell over the 
Fermi energy for two thinnest rings with a) L = 4 (96 
atoms) and b) L = 8 (192 atoms). For smaller ring the 
ground state has total spin 5* = 1 with a very small gap 
to the first excited state with S* = [40|. The lowest 
states with larger total spin have higher energies. For 
192 atoms ring the total spin of the ground state is max- 
imal (5* = 3). The lowest levels with different total spin 
have slightly higher energies. This can be understood in 
a following way. The splitting between levels is large for 
smaller structures, which is seen in Fig. [5l For ring with 
L = 4 (96 atoms) this value (0.17 eV) is comparable with 
electronic interaction terms, e. g. 0.34 eV for two elec- 
trons occupying the lowest state. For ring with L = 8 
(192 atoms) the electron-electron interaction terms are 
around 0.23 eV (for interaction between two particles on 
the first state) , which is much larger in comparison with 
single particle energy difference (0.015 eV). From this, 
we clearly see that for ring with L = 4 it is energetically 
favorable to occupy low energy states by electrons with 
opposite spins. For ring with L = 8 all states have simi- 
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FIG. 6. (Color online) The low energy spectra for the differ- 
ent total spin S of half filled first shell over the Fermi energy 
for two thinnest rings with (a) L = 4 (96 atoms) and (b) 
L = 8 (192 atoms). 



lar energies and due to exchange interactions the lowest 
energy state is maximally spin polarized. 

The behavior of magnetic properties of the ground 
state for half filled shell as a function of size is shown 
in Fig. [71 In this case, ground state properties can 
be explained as a result of the competition between oc- 
cupation of levels with smallest single particle energies 
which favors opposite spin configurations, and parallel 
spin configurations for which exchange interactions are 
maximized. For rings with L > 5 the ground state is 
maximally spin polarized. Here, the splitting between 
levels is relatively small and the ground state is deter- 
mined by electronic interactions. Moreover, this split- 
ting decreases with increasing size and this is seen in the 
spin gap behavior (Fig. [7|). The largest spin gap is ob- 
served for ring with L = 6 and decreases with increasing 
L. For small rings the situation is more complicated. 
Here, the contributions from single particle energies and 
interactions are comparable. As a consequence, we ob- 
serve ground states with alternating total spin S = 1 and 
5* = 0. For sufficiently large rings, L > 5, we observe sta- 
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FIG. 7. (Color online) Lower: Total spin of the ground and 
first excited states for the half filling of the first shell in the 
thinnest ring structures with different sizes. Upper: Corre- 
sponding energy spin gap between ground and first excited 
states. 



bilization of the spin phase diagram. This is connected 
to changes of the energy differences between levels in a 
shell - above a critical size these values are so small that 
they don't play a role anymore. 

In Fig. [8] we show the phase diagram for a ring with 
L = 8 (192 atoms). Near the half filling the ground 
state is maximally spin polarized which is related to the 
dominant contribution from the short-ranged exchange 
interaction terms, and charge density is symmetrically 
distributed in the entire ring (see Fig. 3]). Adding or re- 
moving electrons causes irregularities in the density dis- 
tribution, and correlation effects start becoming impor- 
tant. This results in an alternating spin between maximal 
polarization (e.g. 3, 4, 9 extra electrons) and complete 
depolarization (e.g. 2, 8, 10 extra electrons) of the sys- 
tem. 



VI. CONCLUSIONS 

We presented here a theory of graphene quantum rings 
designed to produce degenerate shells of single particle 
states close to the Fermi level. By combining tight- 
binding and configuration interaction methods, we an- 
alyzed magnetic properties and electronic correlations in 
such structures as a function of size and number of added 
electrons. For the half filling of the degenerate shell in 
sufficiently large ring, maximal polarization of the ground 
state is predicted. Away from the half filling the cor- 
relation effects appear and the ground state total spin 
alternates between maximal polarization and complete 
depolarization. 
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FIG. 8. (Color online) The spin phase diagram for electrons 
occupying first shell of the ring structure with L = 8 (192 
atoms) . 
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APPENDIX A: CALCULATIONS OF COULOMB 

MATRIX ELEMENTS BETWEEN ELECTRONS 

ON SITES IN GRAPHENE HONEYCOMB 

LATTICE. 

The Coulomb interaction term from eq. (J2j), we can 
write as 



^ = n Yl (^P\^ \ df)al^al^,ada'afa 



s,P,dJ, 
a, a' 



2 ^ 

s,P,dJ, 



J2 {ij I V I kl)AlA^^AiA{clc]^,Ck.>ci, 

i,j,k,l, 



"^^^ A^drr and A^ are coeffi- 



where we substitute asa — A^ia^H'-^cr 

cients in transformation from one basis to another. Here 

the Coulomb matrix elements are defined as 

dridr2ilJi (ri)V^| (r2) x 



i^J\V\kl) = — 



1 



r2 - ri 



T^k{r2)^i (ri) 



where hz is an effective dielectric constant and ipi (ri) is 
Slater tTz orbital on a site i of electron 1, given by a 
function 



^(^^) = (|^ 



zexp 



-Cri 



with ^ = 3.14. Below we show selected Coulomb matrix 
elements for /^ = 1. Numbers 1 and 2 and 3 indicate 



{ij 1 V 1 kl) 


E[eV] 


(11 IF 111) 


16.522 


(12 1 y 1 21) 


8.640 


(13 1 y 1 31) 


5.333 


(11 1 y 1 12) 


3.157 


(12 1 y 1 31) 


1.735 


(12 1 y 1 12) 


0.873 


(11 1 y 1 22) 


0.873 


(22 1 V 1 13) 


0.606 



TABLE L Selected coulomb matrix elements between elec- 
trons on sites in graphene honeycomb lattice.. Numbers 1, 2 
and 3 indicate electron on-site, on nearest neighbor site and on 
next nearest neighbor site of hexagonal lattice, respectively. 



electron on-site and on nearest neighbor site and on next 
nearest neighbor site of hexagonal lattice, respectively: 
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